PROCEEDINGS 

0F SCIENCE 



PS 



KEK Preprint 2009-28 
CHIBA-EP-180 



Topological configurations of Yang-Mills field 
responsible for magnetic-monopole loops as quark 



dv confiner 

o 
o 



> 
O 



Akihiro Shibata 

Computing Research Center, High Energy Accelerator Research Organization (KEK) & 
Graduate University for Advanced Studies ( Sokendai), Tsukuba 305-0801, Japan 
E-mail: akihiro . shibata@kek . jp 

Kei-lchi Kondo 

Department of Physics, Graduate School of Science, Chiba University, Chiba 263-8522, Japan 



E-mail: kondok@ faculty . chiba-u . jp 

Seikou Kato 

Takamatsu National College of Technology, Takamatsu 761-8058, Japan 
i-C ■ E-mail: kato@takamatsu-nct.ac.jp 



Shoichi Ito 

Nagano National College of Technology, 716 Tokuma, Nagano 381-8550, Japan 
CO , E-mail: 



CO 

Toru Shinohara 



shoichi@ei . nagano-nct . ac . jp 



Graduate School of Science, Chiba University, Chiba 263-8522, Japan 
E-mail: 



Shinohara@graduate . chiba-u . jp 



Nobuyui Fukui 

Graduate School of Science, Chiba University, Chiba 263-8522, Japan 



E-mail: n . f ukui@graduate . chiba-u . jp 



We have given a new description of the lattice Yang-Mills theory a la Cho-Faddeev-Niemi-Shabanov, 
which has enabled us to confirm in a gauge-independent manner "Abelian"-dominance and magnetic- 
monopole dominance in the Wilson loop average, yielding a gauge-independent dual superconductor pic- 
ture for quark confinement. In particular, we have given a new procedure (called reduction) for obtaining 
a gauge-independent magnetic monopole from a given Yang-Mills field. In this talk, we demonstrate 
how some of known topological configurations in the SU(2) Yang-Mills theory such as merons and in- 
stantons generate closed loops of magnetic-monopole current as the quark confiner, both of which are 
characterized by the gauge-invariant topological index, topological charge (density) and magnetic charge 
(density), respectively. We also try to detect which type of topological configurations exist in the lattice 
data involving magnetic-monopole loops generated by Monte Carlo simulation. Here we apply a new ge- 
ometrical algorithm based on "computational homology" to discriminating each closed loop from clusters 
of magnetic-monopole current, since the magnetic-monopole current on a lattice is integer valued. 
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1. Introduction 

It is interesting to study the confinement mechanism of QCD. The dual superconductivity picture is a 
promising mechanism for quark confinement. It is known that the string tension of the Abelian part and 
the monopole part in Yang-Mills (YM) fields reproduce the original one, so called Abelian dominance and 
monopole dominance in the string tension. There is another approach that the center vortex can explain the 
string tension. However, the dominances of these objects has been observed only in the special gauges such 
as the maximal Abelian (MA) gauge or the maximal center gauge, but not in other gauges. 

We have given a new description of the Yang-Mills (YM) theory on a lattice, which is expected to 
give an efficient framework to explain quark confinement based on the dual superconductivity picture in the 
gauge independent manner. In the case of the SU (2) Yang-Mills theory, it is constructed as a lattice version 
of the Cho-Faddeev-Niemi-Shabanov (CFNS) decomposition!]]]] in a continuum theory [fJ][^][Q]. By per- 
forming numerical simulations we have demonstrated that the gauge-invariant magnetic monopoles can be 
constructed and they reproduces the string tension [g], and shown the infrared "Abelian" dominance in the 
propagator, i.e., the extracted "Abelian" propagator is dominant in the infrared region^]. We have extended 
the framework to the SU (N) YM theory In fact, we have demonstrated the numerical simulation in the 
SU (3) case, which has two options corresponding to its stability subgroups. One is the maximal option 
with the stability subgroup H = t/(l) X £/(l), which is the gauge-independent reformulation of the Abelian 
projection represented by the conventional MA gauge. We have demonstrated the gauge independent study 
of "Abelian" dominance in the conference lattice2007 [JSj] . The other is the minimal one with the stability 
group H = U(2), which is a new type of formulation and derives a non- Abelian magnetic monopole 
We have demonstrated the non- Abelian magnetic monopole dominance in the conference lattice2008 

In this talk, we apply this method to investigate the role of magnetic monopole for the confinement in 
the gauge independent manner. In what follows, we restrict to the SU (2) case. We summarize the result 
of the lattice formulation of the CFNS decomposition, which can extract the dominant degrees of freedom 
that are relevant to quark confinement in the Wilson criterion in such a way that they reproduce almost 
all the string tension in the linear inter-quark potential. Combined with non-Abelian Stokes's theorem, 
the magnetic monopoles can be extracted from the decomposed relevant part of YM field in the gauge 
invariant manner. The monopole dominance suggests that the magnetic monopole plays a central role in 
quark confinement. So, it is interesting to investigate the magnetic monopole as a quark confiner. We 
demonstrate the analysis for some known topological configurations as well as lattice data, where we can 
extract magnetic monopole loops directly from YM configuration, which are the gauge invariant physical 
objects. 

2. New variables on a lattice 



We summarize the new description of the SU (2)YM theory on a lattice [|3|][]4|] (See also [Jllj]): The YM 
field Aj/ jj is represented as a link variable 

/ rx+fle \ 

U x ,ii = exp ( -ig J dx^Ap (x) J = exp (-jgeAy iM ) , (2. 1) 

which is supposed to be decomposed into the product of new variables, X x u, V X jU G SU (2), 

A',,,1',,, (2.2a) 

V x ^ = exp (-igeVy ijU ) , (2.2b) 

X Xtfl = exp (~ige\ x ^) , (2.2c) 
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such that the decomposed variables are transformed by a full gauge transformation Q. x £ SU (2): 



(2.3a) 
(2.3b) 



where V x ^ is defined on a link (jc,x+£fl) like U xfl , and X x ^ on a site. So, the corresponding (Lie algebra- 
valued) gauge filed is evaluated at the midpoint, x' = x+£jj./2, for U x ^ and V x ,^, and at the sitex for X x ,^.To 
define the decomposition, a color field, n x = n x a k /2 is introduced as a site variable, where o k is the Pauli 
matrix and n x (k = 1,2,3) is a unit vector. The color field is transformed adjointly by an independent gauge 
transformation & x G SU (2) as — > & x ii x &l. The decomposition is determined by solving the defining 
equations: 



D e ^[V]n x :=n x V Xjfl 
tr(n x X c# )=0, 

and the solution is obtained in terms of U Xi n and 



V X ,H n X +H 



0, 



1 



f/t 



(2.4a) 
(2.4b) 



(2.5a) 
(2.5b) 



The reduction condition must be introduced in order that the theory written in terms of new variables is 
equipollent to the original YM theory, i.e., the symmetry extended by introducing the color field, SU (2)& x 
[SU (2) /U(l)] Q , must be reduced to the same symmetry as the original YM theory, SU (2)n = ©. The reduc- 
tion condition is given by the minimization of the functional, which is invariant under the gauge transfor- 
mation, SU (2)n = 0. Here, we use the following functional: 

(2.6) 



F Red . = E^tr ( (D^[U]n x ) (D;[U]n x y\ /tr(l) . 



Note that the stationary condition with respect to the color field, dF^ s ±/dn k x = 0, gives the differential form 
corresponding to the continuum theory, and it determines the color field for a given YM field. The algorithm 



solving the reduction condition is given by Ref.[ 11] 



3. Gauge independent magnetic monopole 

The Wilson loop operator on a lattice is given by the path ordered product of link variables, 



W C [U] := tr [pU<xs+h>eC U w] M 1 ) = tt Pex P 



Ml), 



(3.1) 



where the path C is defined along the relevant links, and we have used the definition of the link variable, 
in the 2nd equality. Following the paper the Wilson loop operator in the fundamental representa- 
tion in the continuum theory is rewritten into 



W C [A] :=tr 



Mi) 



Pexp (—ig^dx^A^x 
d^lcexpi-ig^dx^W^ 
fdn[S]zexp\-ig [ dS^ v ^ v [\] 

J I Jl.:d2.=C 



(3.2) 
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where (jc) is given by 



V m (jc) =tr(A M (x)n(x))n(x) + - [^n(x),n(x)] . (3.3) 

^8 



Following the paper [^], it is turn out that the defining equation eq(2.4) gives the decomposition which 



reproduces "Abelian" (V) dominance for the Wilson loop operator on a lattice even with a finite lattice 
spacing £ ; 

Wcp] = consign W C [V] := tr [PU<. w >ec V ^] M 1 ) ^ 

and we can identify the Vy ifl in eq( [2.2b[ ) withV jU (x / ) at the midpoint of the link. The field strength of V is 
given by 

^,« y *+«,j3 V,Vj3,« y i = ex P H# £tr («^.Mv t V K) n y ) . (3.5) 



By using the Hodge decomposition, eq(p.2D is further rewritten into 



W C [A] = J dn[Z]zexp{-ig(K,Zz)-ig(J,Nz)}, (3.6) 

where we have used K := 8*F, J := 8F for F being the field strength 2-form J^ v , and := 5*SxA _1 and 
Ny. '.= 8&zA . with the four-dimensional Laplacian, A = d8 + Here 5x is the vorticity tensor defined 
by s£ v = J z dS^ v (X(o))8(x-X(o)) on the surface L:dL = C supported by the Wilson loop C. Therefore, 
the monopole dominance for the lattice Wilson loop is given by 

(W C [U\) * (W C [V]) = (exp (-igeZ^^)) > ( 3 - 7a ) 
where the lattice magnetic monopole current K x ^ and S v ^ are given by 

E x ^=Y, s ^ap 7 da^ 1 {x-s)S J p7 {s + ^l), d a S J ap =J p GC (3.8) 

k x ^ = 2xk x ,n = ^s^apd^p [V] , ® x>a p [V] := argtr ((1 + n^aVv+a^W^) , (3.9) 

where d% denotes the forward difference (lattice derivative) in the X direction: ed^f(x) :=f(x + eX)—f(x). 
It should be noticed that the magnetic monopole current k x ^ is gauge invariant, as can be seen from the 
transformation law of the new variables. 

4. Analysis of the magnetic monopoles as quark confiners 

We have demonstrated the "Abelian" dominance and the magnetic-monopole dominance in the gauge 



invariant wayjl 1|]. The static potential calculated from the decomposed variable V is reproduced by one 
from the Yang-Mills field, i.e., the string tension obtained by V field 94% of one from YM field. The string 
tension calculated only from magnetic monopoles reproduces 93% of it. This shows that the proposed 
new description of YM theory extract the degrees of freedom relevant for the quark confinement from YM 
field, and the magnetic monopole which is the decomposed part of V x ilx , plays a central role for the quark 
confinement. So we investigate the magnetic monopoles as quark confiners. 

The magnetic monopole current is a integer- valued, k X n € {0, ±1,±2}, and can be identified with the 
link variables on the dual lattice (x — ejl,x) with x = {x\ + e/2, X2 + e/2, xj + e/2, X4 + e/2). So non-zero 
monopole currents can be identified with directional edges which construct a directional graph. The current 
conservation, 
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Figure 1: (Left panel) The plot of a magnetic-monopole loop generated by a pair of (smeared) merons in 4- 
dimensional Euclidean space. The 3-dimensinal plot is obtained by projecting the 4-dimensional dual lattice space to 
the 3-dimensional one, i.e., (x,y,z,t) — > (y,z,t). The positions of two meron sources are described by solid boxes, and 
the monopole loop by red solid line. (Right panel) The plot of the topological charge density for z — t plane (slice of 
x — y = 0). Two peaks of the topological charge density are located at the positions of two merons. 



implies that there is no source and no sink. When we identify a monopole current of charge k x ^ = ±2 with 
a double edges of a single charge, each vertex (site of dual lattice) has the same number of incoming and 
outgoing edges. Therefore, the monopole currents construct the loops. Here, some loops are connected 
each other at dual lattice sites, or share links of the dual lattice carrying the double monopole charge. Thus 
the analysis of the monopole configuration is converted to a geometrical problem. In the magnetic part of 



the Wilson loop operator eq(3.6) or eq( 3.7a ), the term (K,Ez) represents the 4-dimentional version of the 
Gauss's linking number formula, that is, the linking number between a single loop K and a two-dimensional 
surface £ whose boundary is the Wilson loop C. The conjecture on translation the contribution of center 
vortex and of magnetic monopole loops to the Wilson loops in Ref. |Q] can be checked by using the lattice 
data. 

Now, let us apply this method to topological configurations of Yang-Mills field, which has been con- 
sidered to play the central role for confinement. The topological solution can be translated to the lattice 
variable by using the definition of the link variable: 

= Pexp (-igj* + Wa m (x)) - ^nLi ex P (-^ A m(*+ (« " V*)^)) , 

where the integral along the link (x,x + jle) is calculated by the path ordered product of the exponenti- 
ated gauge potential. Here, we demonstrate the two-merons case. Figure [j] shows the detected magnetic 



monopole. The numerical analysis on the lattice reproduces the result of analytical study in Ref.[|12|], even 
though the analysis is done by using the finite volume lattice. It is interesting to investigate other topo- 
logical configurations such as two-instantons and calorons. We find a monopole loop for the JNR type of 
two-instantons solution, while no monopole loop for one-instanton and two-instantons of 't Hooft type. The 
detail analysis of two-instantons solutions will appear in the subsequent paper [|l3|]. 

Then, we return to the analysis of lattice data. It is very hard to manipulate monopole configurations 
directly, since they contain more than 35000 non-zero monopole currents (see the left panel of Figure ||). 
Therefore, we introduce the algebraic algorithm for topology. The CHomP homology software, provided 



by the computational homology project[14], computes topological invariants called the Betti numbers of 
a collection and their generators in the algebraic way. The Betti numbers are part of the information con- 
tained in the homology groups of a topological space, which intuitively measure the number of connected 
components, the number of holes, and the number of enclosed cavities in low dimensions. In our case, the 
generators of the dimension-one homology group correspond to magnetic monopole loops. 
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Figure 2: (Left panel) The analysis of magnetic monoples for 400 configulations in 24 4 lattice with the parameter 
P = 2.4. The blue line shows the number of the non-zero magnetic monoole currents (right virtical axes). The red 
line an dgreen line show the number of clusters of connected loops (the Betti number of dimension zero) and the 
number of loops (the Betti number of dimension one) for each configuration, respectively (see left virtical axes). 
(Right panel) The 3-dimensional polt of detected magbetc monopole loops, where the graph in 4dimensional Euclid 
space is projected to the 3-dimensional space, i.e., (x,y,z,t) — > (x,y,z). 




Figure 3: (Left panel) The histgram for the length of monopole loops, i.e., for each configuration the number of 
monopole loops with length n is counted, and is polted as a truss of poles with the same length. (Right panel) The 
histogram of the cluster size. 

We apply the method to the lattice data of 24 4 lattice with periodic boundary condition whose config- 
urations are generated by using the standard Wilson action with the parameter j8 = 2.4. The left panel of 
Figure ^ shows data of detected magnetic monopole currents for 400 configurations. The blue line shows 
the number of non-zero charge currents, i.e., non-zero magnetic monopole currents share about 3% of links. 
The red line and green line show the number of clusters and the number of loops for each configurations, 
respectively. The right panel of Figure |2| shows an example of detected magnetic monopoles, which are 
plotted in the 3-dimensional space projected from the 4-dimendional Euclidian space. Figure [3] shows de- 
tail of the magnetic monopole configurations. Then, we are ready to investigate the magnetic monopole 
contribution to the static potential by using extracted monopole loops. 

5. Summary and discussions 

We have given a new description of Yang-Mills theory on a lattice, which gives the decomposition of 
the Yang-Mills field for extracting the relevant degrees of freedom for quark confinement in the gauge in- 
dependent manner based on the dual superconductivity picture. The extracted magnetic monopoles explain 
the string tension and they should play the central role for quark confinement. It is interesting to investigate 
the magnetic monopole as a quark confiner. The implication of the quark confinement by the configurations 
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topological charge can be convert to the geometrical relations of the topological charge density and distribu- 
tion of the magnetic monopole loops, since the magnetic monopole and topological charge density are given 
as a gauge invariant object (physical object). We have applied the method to the topological configurations 
and the lattice data, and extracted the magnetic-monopole loops directly. Using these magnetic-monopole 
configurations, we investigate the relations among the topological charge density, magnetic monopoles and 
static potential. It is also interesting to investigate the implication of the magnetic monopole loops for the 
confinement and deconfinement phase transition. 
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